A Note on Overshoot Estimation in Pole Placements 
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In this note we show that for a given controllable pair (A, B) and any A > 0, a gain matrix K can be chosen so that the 
transition matrix eS A+BK ^* of the system x = (A + BK)x decays at the exponential rate e~ xt and the overshoot of the 
transition matrix can be bounded by MX L for some constants M and L that are independent of A. As a consequence, for any 



h > 0, a gain matrix K can be chosen so that the magnitude of the transition matrix e 



{A+BK)t 



can be reduced by | (or by 



any given portion) over [0,h]. An interesting application of the result is in the stabilization of switched linear systems with 
any given switching rate (see [1]). 
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1 Introduction 

Consider a linear system 



x = Ax + Bu, 



(1) 



where x(-) takes values in R™, it(-) takes values in K m , 
and where A and B are matrices of appropriate dimen- 
sions. Suppose (A, B) is a controllable pair. It is a well 
known fact that for any A > 0, a gain matrix K can 
be chosen so that the transition matrix of the system 
x = (A+ BK)x decays exponentially at the rate of e _At , 
that is, for some R > 0, 



^(A+BK)t 



< Re 



for some constants M and L independent of A. As a con- 
sequence, one sees that for any h > 0, a gain matrix K 
can be chosen so that the magnitude of the transition 
matrix g( j4 +- B - R ") t can be reduced by ^ (or by any given 
portion) over [0, h\. Note that this is a stronger require- 
ment than merely requiring e(" 4+B - ftr ) t to decay at an ex- 
ponential rate. An interesting application of the result 
is in the stabilization of switched linear systems with a 
given switching rate (see [1] ) . 

The estimate of the overshoots of transition matri- 
ces in the practice of pole assignments has been studied 
widely (see e.g. [5], [9] and [7]). Our main result in this 
note can be considered an enhancement of the Squash- 
ing Lemma (see [7], [6] and [4]) which says the following: 
for any To > 0, 6 > 0, any A > 0, it is possible to find K 
such that 



where and hereafter || • || denotes the operator norm in- 
duced by the Euclidean norm on M. n . To get a faster de- 
cay rate, it is natural to consider a "higher gain" matrix 
K\. However, such a gain matrix in general results in a 
bigger overshoot for the transition matrix e (A+BKi)t_ j n 
this note, we show that in the pole placement practice, 
a gain matrix K can be chosen so that the overshoot of 
the transition matrix e (- 4 + Bi:s ")* can be bounded by MX L 



D (A+BK)t 



< 5e- x{t - To) . 



(2) 



In the current note, we show that K can be chosen so 
that the estimate in (2) can be strengthened to 



M+BK)t <M yL, 
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for some constants M and L which are independent of 
A. Our proof is constructive that shows explicitly how 
M and L are chosen. 



2 Main Result 



tion 



In this section we present our main result. 
Proposition 2.1 Let A e R nxn and B e R nxm be two 
matrices such that the pair (A, B) is controllable. Then 
for any A > 0, there exists a matrix K G flj mx ™ ^ ^at 



e (A+BK)t < MX L e -\t^ Vf > 0; 



(3) 



where L = (n — l)(n + 2)/2 and M > is a constant, 
which is independent of A and can be estimated precisely 
in terms of A, B and n. 

Comparing with the Squashing Lemma obtained in 
[7], Proposition 2.1 has two improvements: (i). In (2), 
the estimate on the transient overshoot is exponentially 
proportional to the decay rate A, which resulted in an 
estimation of the transition matrix in terms of e -Mt- T o) 
instead of e _At . In (3), the estimate on the transient over- 
shoot is proportional to X L instead of e Xr ° as in (2). This 
distinction between the two types of estimations may be 
significant for some possible extensions of our results to 
systems with external inputs, (ii). The value of the con- 
stant M in estimate (3) can be precisely calculated by 
using our constructive proof (see equation (10) in the se- 
quel) . This is certainly a very desirable feature for prac- 
tical purposes. See Example 3.1 for some illustrations. 

Proposition 2.1 was primarily presented and applied 
to a stabilization problem of switched linear systems in 
[2] . It was found later that a recent paper [3] also provides 
a similar result with similar proofs. The difference is that 
[3] only considered the single input case and the upper 
bound MX L in (3) was found to be a polynomial p(X) in 
[3] without an explicit expression. Hence, our result has 
obvious merits in control design. 

Proof of Proposition 2.1. First we consider a linear 
system (^4, b) of a single input. Without loss of general- 
ity, we assume that (A, b) is in the Brunovsky canonical 
form: 



A = 



/o 1 ••• \ 
1 ••• 

\ai a 2 a 3 ■■ ■ a n J 



M 
o 

w 



Let Ai,...,A„ be n distinct, negative real numbers. 



There exists some k e 



such that the charac- 



teristic equation of the closed-loop system A + bk is 
p(X) = (A - Ai)(A - A 2 ) • • • (A - A„). Note that the 
closed-loop system is given by 

x\ = x 2 , x 2 — x 3 , x n -i = x n , 

x n = Pixi + (3 2 x 2 H (3 n x n 

for some (3 2 , . . . , (3 n £ R. Hence, x\ satisfies the equa- 



(n) 



PlXl + f3 2 Xi H h P n x{ 



(n-1) 



(4) 



whose characteristic equation is the same as p( A) . Hence, 
the general solution of (4) is 



xi(t) = cie Alt + c 2 e 



,X n t 



where c\ , c 2 , . . . , c n are constants. From the equations 
x 2 = ±i,x 3 = x 2 ,...,x n — x n -\, we have x(t) = 
A e Dt c, where 



A n = 



/ 1 1 

Ai A 2 



1 \ 

A n 



D - 



\r 1 : 

Ai o 

o A 2 



\n-i \n-l / 







\o 

and where c 



XnJ 



^ci c 2 ■ ■ ■ c^j ■ Now, observe that 

x(0) = Aoc, that is, c = Aq 1 x(0) (note that Ao is an 
invertible Vandcrmondc matrix). Comparing this with 
the transition matrix of the system, one sees that 



e (A+bk)t =Aoe Dt A -l_ 



(5) 



Let A max = max{|Ai| , . . . , |A„|}. Without loss of gen- 
erality, assume that A max > 1. To get an estimate on 
||Ao|| and ||Aq ||, we need the following simple fact: for 
an n x n matrix C, let c max = maxi<i.j<„ |. It is not 
hard to see that 

l|C||< 

Hence, we have 

||A ||<nA^. (6) 

To get an estimate on A \ first note that 
1 



An" 1 



■ adj A , 



(7) 



detA 

where adj Ao denotes the adjoint matrix of Ao, and that 
detAo = Y[{Xj - A»). 



j>i 



Hence, if we choose Ai, . . . A„ in such a way that Ai + i < 
Aj - 1 with Ai < 0, we get |det A | > 1. 



2 



Taking the structure of adjAo into account, it is easy 
to see that for C = adjA , 

c <fn-lVA i+2+-+(n-i) 

"-max _^ V / • ' v max 

n(n-l) 

= (n-l)!A max 2 . 



Hence, by (7), we have 

IIAq 1 !) < ||adjA || <n(n -1)!A D 
Consequently, (6) and (8) yield that 



'■■»-!) 



(8) 



||Aoe Dt A^||<nA^ || e Dt || n(n - 1)!A^" 1)/2 
^nnlA^" 1 )^ 2 )/^-^*, 

where A min = min{|Ai| , . . . , |A„|}. 

Suppose for some p > 1, A max < pA m i n - Then, it fol- 
lows that 



\A e Dt A^\\ < MA^ 1)( " +2)/2 e- Am '°*, 



(9) 



where 

M = nn\p i - n - 1 ^ n+2 ^ 2 . (10) 
In summary, we need the following conditions on the Ai 's: 

• Ai, A 2 , • • • , A„ are distinct, real, and negative; 

• Aj + i < Ai — 1 for 1 < i < n—1, and hence, A max = A„|, 

A m in = | Al |; 

• |A„| < p | Ai|, for some constant p > 1. 

Obviously, for any given A > 0, it is easy to choose 
Ai, • • • , A n to satisfy all the above conditions together 
with the condition that Ai < —A. For example, one can 
choose Ai < min{ — 1,— A}, and let Ai+i = Aj — 1 for 
1 < i < n — 1. Since |A„| = |Ai — (n — 1)| < n \X\\, we 
see that p can be set as p — n. 

With such choices of Ai, A 2 , . . . , A„, we see from (5) 
and (9) that the desired result hold. 

Now we consider the case when (A, b) is not in the 
Brunovsky canonical form. In this case, find an invertiblc 
T e W nxn such that (T^AT, T _1 6) is in the Brunovsky 
canonical form. 

For any given A > 0, the above proof has shown that 
for A 1 = T^AT, h = T^b, one can find fc £ R lx ™ 
such that 

e (A 1 +b 1 k )t < MX L e~ Xt , 

where M is given by (10) for some chosen p, and L = 
(n - l)(n + l)/2. Clearly, with k = fcoT" 1 , one has 



,{A+bk)t _ 



= T(e 



(Ai+6ife )t 



)T" : < MiA L e 



L-\t 



where M 1 = M \\T\\ \\T~ l \\. 

Finally, we consider the multi-input system 

x = Ax + Bu, 



(11) 



(12) 



where i £ R" X ",B e R" xm . Suppose that the system 
is controllable. By Heymann's Lemma (c.f., e.g., page 
187 of [8]), one sees that for any v e R m such that 
b := Bv ^ 0, there exists some K G R mx ™ suc h that 
(A + BKo, b) is itself controllable. Hence, the conclusion 
of single-input case that has just been proved above is 
applicable to the controllable pair (A + BK 0} b), and 
one then sees that there exists some k € R lx ™ such that 
|| e (A+BK +bfc)t|| < MX L e~ xt for all t > 0. Hence, with 
K = K + vk, it holds that 



e (A+BK)t <MX L e -\t vt>()> 



(13) 



This completes the proof. □ 

Remark 2.2 In the above proof, we have used the fact 
that for a single input system (A, b) which is controllable, 
when it is not in the Brunovsky canonical form, one can 
find an invertible matrix T such that (T^AT.T^b) is 
in the canonical form. To be more precise, the matrix T 
can be chosen as (sec e.g., [8]): 



b Ab 



/ a„_i • • • ai l\ 

: : 1 

ai 1 ••• 

\ 1 ••• 0/ 



where oi, . . . , a„_i are as in the characteristic polyno- 
mial of A given by 



det(sI-A) =s n +a 1 s n - 1 + 



+ fl n - 2 s + a n -i. 



From this one can find an estimate of ||T|| and ||T 1 ||, 
which in turn will lead to an estimate of Mi in (11). 

3 An Example 

The design technique is demonstrated in the following 
example. 

Example 3.1 Consider the following controllable linear 
system: 



A 1 



(10 l\ 

1 -1 

V 21 o/ 



B x 





w 



With the help of MATLAB, we first calculate the trans- 
fer matrix 

/ l\ 



Ti 



-1 -1 

V-i o v 
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With the transfer matrix T\ , one has 



/ 1 o\ 

1 
V-l 2) 





W 



Calculation shows that ||Ti|| = 1.80193754431757 and 
HTf 1 !! = 2.24697960199992. Taking p = n(= 3), we 
have 



L _ (n-l)(n + 2) _ 5 



(14) 



M = ||Ti||||T 1 - 1 ||nn!n( n - 1 )( n - 2 )/ 2 « 218.642. (15) 

Suppose for some design purpose, a decay constant 
A = 49.894 is given. Choosing Ai = —A, A 2 = Ai — 1, 
A 3 = A 2 — 1, the feedback K\ can be easily calculated 
(under the normal form) as 



-151.681 -7769.474 -131773.562 



)• 



Back to the original coordinate frame, we have 



124155.769 7769.474 -7617.793 



■)■ 



With such a choice of Ki, we get the desired decay esti- 
mate 



,(A+BK)t 



< M\ L e' xt Vf > 0, 



for the given decay constant A = 49.894 with L and M 
given as in (14)-(15). □ 

4 Conclusion 

In this note we show that if (A, B) is controllable, 
then for any A > 0, a gain matrix K can be chosen such 
that the transition matrix g( A + BK ) t decays at the ex- 
ponential rate e~ xt and the overshoot of e ( A + BK ) 1 can 
be bounded by M\ L for some constants M and L that 
are independent of the decay constant A. The result pro- 
vides a convenient tool for control design, particularly 
for switched systems, see [1]. 
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Erratum 

There is a mild flaw in the statement of Proposition 
2.1 in the above paper (cf. [1]). We restate it as follows. 

Proposition Let A e R nxn and B e R" xm be two 
matrices such that the pair (A, B) is controllable. Then 
for any A > 1, there exists a matrix K 6 R mx ™ suc h that 



e (A+BK)t < MX L e 



-\t 



Vi > 0, 



(16) 



where L = (n — l)(n + 2)/2 and M > is a constant, 
which is independent of A and can be estimated precisely 
in terms of A, B and n. 

The proof of Proposition 2.1 in [1] is only valid for 
the case when A > 1 (instead of the original version of 
A > 0), because the eigenvalues Ai, . . . , A„ were chosen 
to satisfy Ai < —1, and Afe < Ai for k > 1. For more de- 
tails, we refer the reader to the discussions that followed 
formula (10) in [1]. 

A main motivation of the work in [1] was for us to 
develop the results in [2]. As in most applications of 
overshoot estimation for pole placements, the parameter 
A in [2] was chosen as a number of large value. Hence, 
the correction does not affect our results in [2] . 

Acknowledgment. The authors would like to thank 
Prof. Elena De Santis for pointing out the error. 
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